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ABSTRACT 



The origin of SO galaxies is discussed in the framework of early mergers in a Cold Dark Matter 
cosmology, and in a scenario where SOs are assumed to be former spirals stripped of gas. From 
an analysis of 127 early-type disk galaxies (SO-Sa), we find a clear correlation between the scale 
parameters of the bulge (r e ff) and the disk (/ir), a correlation which is difficult to explain if these 
galaxies were formed in mergers of disk galaxies. However, the stripping hypothesis, including 
quiescent star formation, is not sufficient to explain the origin of SOs either, because it is not 
compatible with our finding that SOs have a significantly smaller fraction of bars (46±6 %) than 
their assumed progenitors, SO/a galaxies (93±5 %) or spirals (64-69 %). Our conclusion is that 
even if a large majority of SOs were descendants of spiral galaxies, bars and ovals must play an 
important role in their evolution. The smaller fraction particularly of strong bars in SO galaxies 
is compensated by a larger fraction of ovals/lenses (97±2 % compared to 82-83 % in spirals), 
many of which might be weakened bars. We also found massive disk-like bulges in nine of the SO 
galaxies, bulges which might have formed at an early gas-rich stage of galaxy evolution. 

Subject headings: galaxies: evolution — galaxies: bulges — galaxies: structure 



1. Introduction 

The formation of SO galaxies is generally dis- 
cussed in the framework of mergers of disk galax- 
ies within a Lambda Cold Dark Matter (ACDM) 
cosmology. Alternatively, they can be viewed as 
former spirals where star formation has ceased and 
gas lost by stripping mechanisms. However, less 
attention has been paid to the role of bars in the 
evolution of these galaxies. Bars are known to be 
efficient drivers of gas towards the central regions 
of galaxies, and in th e presence of nuclear bars 
( Shlosman et al. 1989) or nuclear spirals, a cen- 
tral starburst might occur. Numerical simulations 



(lAthanassoulal 12003 



Martinez- Valpuesta et al 
20061 ) also predict that bars evolve due to angular 
momentum transfer between the bar and a massive 
or centrally concentrated halo, leading to more 
prominent bars. Indeed, there is observational 
evidence that bars in early-type disk galaxies 
might be more evolved t han bars in spiral galax- 
ies. They are longer (lElmegreen &: Elmegreen 



1985; Laurikainen et al 



20071) . 



more massive. 



and have more freq uently double-peaked Fo urier 



amplitude profiles ( Laurikainen et al. 2007t) and 



ansae-type morphologies dLaurikainen et al. 2007 



Martinez- Valpuesta et al.1 120071 ). In the above 
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simulation models these features can be inter- 
preted as indices of evolved bars. 

Lenses are features commonly observed in SO 
galaxies, but are sometimes also seen in early-type 
spirals. They are defined as components with a 
shallow or constant surface brightness profile, and 
a sharp outer edge. Lenses are a fu ndamental par t 
of the original classification of SOs (|Sandagdll96lh , 
although they were not recognized initially with 
their own type symbol. Ovals are global devi- 
ations of the disk from the axisymmetric shape. 
In distinction to bars they have lower ellipticities 
and generally lack higher order Fourier terms than 
m=2. 

We emphasize the importance of ovals and 
lense s for unders t andin g bar-induced galaxy evolu- 
tion. Kormendv ( 1979h showed that 54% of barred 
galaxies of types SBO-SBa have lenses. The ac- 
tual frequency of lenses in non-barred SOs was not 
determined at that time, although st rong lenses 
were known to exist in such ga laxies (jKormendv 
1984tlSandage fc Brucat"olll979h . Our recent stud- 
ies (Laurikainen et al. 2005, 2007) have confirmed 
that lenses are co mmon in non - barre d SOs. It was 
first suggested by Kormendvl ( 19791 ) that lenses 
might be destroyed or weakened bars. If bars 
evolve in the Hubble sequence, and lenses are in- 
deed weakened bars, at some stage one would ex- 
pect the fraction of lenses to exceed the fraction 
of bars. It is also challenging to explain the origin 
of the multiple bars, ovals and lenses seen in sin- 
gle SOs and for which no explanation has yet been 
given in the current paradigm of galaxy formation. 

In this Letter the origin of SO galaxies is dis- 
cussed, based on an analysis of 127 early-type 
disk galaxies (Laurikainen et al. 2005, 2006): 
N(S0)=82, N(S0/a)=18, N(Sa-Sab)=27. Al- 
though by far most SOs might be descendants 
of spiral galaxies, we present evidence that bars, 
ovals and lenses have played an important role in 
their structure formation and evolution. Support 
for the stripped spiral hypothesis is provided by 
a clear correlation between bulge and disk scale 
parameters (Section 4). The role of bars in struc- 
ture formation is evidenced by the large number 
of ovals+lenses (interpreted as weakened bars) in 
SOs, and by the massive disk-like bulges found 
in nine of the galaxies (Sections 3 and 5). The 
studied galaxies are part of a magnitude-limited 
Near-IR SO galaxy Survey (NIRS0S), which has 



the following selection criteria: Bt < 12.5, in- 
clination < 65°, a nd Hubble type —3 < T < 1 



(|Laurikainen et alJliool iButa et alJl2006t ). More 



than half of the 184 NIRS0S galaxies are currently 
analyzed, which form an adequate sample to study 
the structural components of these galaxies. 

2. Two examples of oval/lens-dominated 
galaxies 

We show 2-dimensional decompositions of the 
surface brightness distribution for two oval/lens- 
dominated galaxies, with the following aims: (1) 
to show how to identify lenses in galaxies, (2) 
to stress the importance of accounting for bright 
ovals/lenses while deriving the parameters of the 
bulge, and, (3) to demonstrate the complexity of 
some SO galaxies that needs to be explained by 
galaxy evolutionary models. Our decompositions 
use a Sersic function for the bulge, allowing for its 
flattening and deviation from elliptical isophotes, 
an exponential function for the disk, and either a 
Sersic or Ferrers function for the bars, ovals and 
lenses. In this study two types of bulges are con- 
sidered: (1) classical elliptical-like bulges (with 
Sersic index n near 4), and (2) disk-like pseudob- 
ulges (with smaller Sersic index), formed mainly 
from the disk material via central star formation. 
However, the vertically thick boxy/peanut struc- 
tures in barred gala xies (often also called bulges) 
(|Athanassoulall2003l ) are considered here as part 
of the bar. Such structures, and nuclear bars and 
rings are not included to the flux of the bulge in 
our decompositions. 

The galaxies and their decompositions are 
shown in Figures 1 and 2: NGC 524 is dominated 
by two almost circular lenses, whereas NGC 5365 
has two oval-shaped components, both ovals with 
an embedded bar. The lenses in NGC 524 are 
directly visible in the image, and showing in the 
surface brightness profile as distinct exponential 
sub-sections. The best fit is obtained by n=2.8 
Sersic bulge and two flat Ferrers functions (Fer- 
rers index=l) for the lenses, implying a fairly 
small bulge-to-total flux ratio B/T=0.28. For 
NGC 5365 the inner bar/oval system is fitted by a 
single Sersic function and the outer one with a sin- 
gle Ferrers function, leading to a fairly exponential 
bulge (n=2.0) and a small B/T=0.17. Counting 
all the flux above the exponential disk as a bulge, 
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would lead to a considerable overestimate of the 
bulge flux (B/T ~ 0.5 in two component fits), in 
accordance with Laurikainen et al. (2005) who 
showed that by omitting the bar/oval in the de- 
composition the B/T-flux ratio is overestimated, 
regardless of whether fD or 2D-decompositions 
are used. It is challenging to explain how this 
kind of multiple ovals/lenses form in SO galaxies 
and how they are related to the evolution of bars. 
In Section 6 they will be discussed in the context 
of cosmologically motivated simulations by Heller 
et al. (2007). 

3. Weak bars inside the lenses 

Although bars in SO galaxies are on average 
fairly prominent, weak bars are detected inside 
the ovals/lenses in eight of the galaxies: NGC 
484, NGC 507, NGC 1161, NGC 1351, NGC 
2768, NGC 2902, NGC 3998 and NGC 7377. All 
these galaxies are classified as non-barred (de Vau- 
couleurs et al. 1991; RC3) and no bar is directly 
visible in the -Kg-band image. However, a weak, 
genuine bar is visible in the residual image after 
subtracting the bulge model obtained from our de- 
composition (see Fig. 3 for NGC 3998). Taking 
into account that the evolution of bars and bulges 
might be coupled, it is interesting to look closer at 
the properties of the bulges of these galaxies. 

We find fairly small Sersic indexes for the bulges 
in these nine galaxies (< n >=2.5), similar to the 
typical values recently found for SO galaxies in de- 
compositions where a multi-component approach 
is used (Laurikainen et al. 2005, 2007; Gadotti 
2008). However, their average bulge-to-total flux 
ratio (< B/T > = 0.44) is higher than average for 
their Hubble type, which for SO-SO/a galaxies ~ 
0.25-0.28 (Laurikainen et al. 2005, 2007; Gadotti 
2008). In Section 6 we will discuss that the weak 
bars, and probably also the massive, fairly expo- 
nential bulges in these galaxies, might be a mani- 
festation of bar-induced secular evolution. 

4. Comparison of the scale parameters of 
the bulge and the disk 

The scale parameters of the bulge and the 
disk are sensitive to the evolutionary processes of 
galaxies during their cosmic history, and therefore 
offer an independent test of the importance of sec- 
ular processes in SOs. A correlation between the 



scale parameters is expected in models where the 
bulges were formed in slow secular processes (re- 
viewed by Kormendy & Kennicutt 2004) , whereas 
no correlation is expected in models where the 
bulges were formed eithe r in a fast dissipative col- 
lapse (lEggen et al J 19621). o r by mergers of large 
( Toomre fe Toomrd Il972t ISpringel fc Hernauistl 
20051 ) or small galaxies (|Abadi et al.ll2003l ). 

The radial scale length of the disk, /ir, and the 
absolute magnitude, M, in the if s -band are calcu- 
lated in the following manner: 

/lR,mtr [kpc] = /lR : obs X D[kpc] / Ci . 

M intT = m ohs - m oxt - c 2 - (5 x log D[Mpc]) - 
25.0, 

where the internal dust correc tion for ft.R, Ci = 
[1.02 - 0. 13 x log (cos i)] , is from l Graham fe Worlev 
(|2008l ). For M, the internal dust correction is 
[0.11 + 0.79 x (1-cos i) %77 ], taken from 



c 2 



Driver et al. (|2008l) . For galaxy distance D we 
use H a — 75 km s" 1 Mpc " 1 , and, Galactic ex- 



tinction, m oxt , is from ISchlegel et al.l (|1998l ). For 
the effective radius of the bulge, r e ff, no inter- 
nal dust correction was made. The subscripts obs 
and intr refer to measured and intrinsic values, 
respectively. We find a clear correlation between 
r c s and /iR^ntr, with the coefficient of correla- 
tion of 0.66, at the significance level of 7E-14 
(Fig. 4). The correlation is independent of the 
applied corrections. The correlation is found to 
be the same for bright (M < -24.5) and faint 
(M > -24.5) galaxies. Thi s cont radicts the re- 
cent result by Bar way et alj ( 20071 ) who found an 
anti-correlation for bright SOs, and a positive cor- 
relation for faint SOs. This different result is most 
probably due to the deeper images and the more 
homogeneous database used in the present study. 
Notice also that the scale parameters in Barway 
et al. were derived using simple bulge-disk de- 
compositions, whereas we use a multi-component 
approach. The correlation we find for SOs is sim- 
i lar to that found pr e viously for sp i ral ga laxies 



(|Courteau et alj Il996t ICarollo et all 120071) . and 



for 14 SB0 galaxies (|Aguerri et al.ll2005l ). 



5. Frequency of bars, ovals and lenses 

Finally, we compare the fractions of galaxies 
with bars and ovals/lenses in different Hubble 
types. If SOs were simply stripped spirals one 
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would expect similar bar fractions in SOs and 
in their spiral progenitors. In Table 1 we use 
RC3 family classes for calculating bar fractions 
in different Hubble type bins: for SO-SO/a galax- 
ies the whole NIRSOS sample of 184 galaxies is 
used, whereas for spirals we use the similarly sized 
Ohio State University Bright Spiral Galaxy Sam- 
ple (OSUBSGS, Eskridge et al. 2000). We find 
that SOs (46±5 %) have bars (SB+SAB) less fre- 
quently than SO/a galaxies (77±9%) or spirals 
(61-70 %). The values for spirals are in agree- 
ment with those found pr e viously by other au- 
thors dKnapen et all \200& lEskridee et all l2000t 



Laurikaincn et al. 2004; Menhdez-Delmestre et al 



2007; iMarinova fc Jogedl2007f ). The bar fractions 
for SO and SO/a galaxies in the sub-sample of 
127 NIRSOS galaxies (38±5 and 76±10 %, respec- 
tively) are nearly the same as for the complete 
NIRSOS sample. In Table 1 the fractions of multi- 
ple bars are calculated in respect of the total num- 
ber of barred galaxies, while all the other values 
are given in respect of the total number galaxies 
within the Hubble type bin. The uncertainties are 
estimated from Ap — y/ (1 —p)p/N, where p de- 
notes the fraction in question in a sample of N 
systems. 

We then use the A2, the maximum m=2 Fourier 
density amplitude in the bar region, normalized to 
m=0, to study three bar strength bins within each 
Hubble type bin. For the OSUBSGS we use the 
values from Laurikainen et al. (2004) , whereas for 
SO and SO/a galaxies they were calculated in this 
study in a similar manner. As a lower limit for 
the barred galaxies we use ^2=0.1. We confirm 
the above result that SO galaxies have a smaller 
fraction of bars than SO/a galaxies or spirals. We 
find that: (1) Sc-Scd spirals have the largest num- 
ber of weak bars (A2 — 0.1-0.3) and the smallest 
number of strong bars ((A2 > 0.6), and that the 
fraction of strong bars increases towards the SO/a 
galaxies. This is in agreement with the previous 
studies showing that the prominence of bars in- 
creases towards the early-type disk galaxies. (2) 
Quite interestingly, although the fraction of strong 
bars increases from late- type spirals towards SO/a 
galaxies, it suddenly drops from 38±9 % to 10±3 
% for SOs. 

We also find that (3) SO galaxies have a larger 
fraction of ovals/lenses than SO/a galaxies galax- 
ies (97±2 v.s 82±9 %), and that (4) SO/a galaxies 



have a larger fraction of multiple bars than earlier 
or later type galaxies (see Laine et al. 2002; Er- 
win & Sparke 2002). The fraction of ovals/lenses 
is found to be the same for barred and for non- 
barred SO-SO/a galaxies (82±4 % v.s 86±6 %, re- 
spectively). Although ovals and lenses might have 
different light distributions, they are not distinct 
enough to be considered separately in our statis- 
tics. 

6. Discussion and conclusions 

In the current paradigm of galaxy formation, 
ACDM, the spheroidal components of galaxies 
were formed through mergers of disk galaxies: dry 
mergers are suggested to lead to the formation of 
elliptical galaxies, whereas mergers of gas poor 
with gas rich galaxies lead to the formation of 
bulges in the disk-dominated galaxies (Khockfar 
& Burkert 2003; Naab et al. 2006). Minor merg- 
ers are actually more common in the Universe 
and they are suggested to form even 55 % of the 
spheroid stars from accreted satellites (Abadi et 
al. 2003). In ACDM the disks form after a major 
merger when hot gas in the halo settles into the 
disk (Kauffmann et al. 1999; Springer & Hern- 
quist 2005). In this picture every dark matter 
halo is expected to possess a substantial pressure 
supported classical bulge with elliptical-like pho- 
tometric properties (Steinmetz & Navarro 2002). 
Although many observations support this scenario 
it is also faced with severe problems, for example 
the bulges, not only in spiral galaxies, but even 
in SO-galaxies are fairly disk-like and have smaller 
bulge-to-total mass flux ratios than predicted by 
cosmological models. Our finding that the scale 
parameters of the bulge (r e s) and the disk (TirJ 
are well correlated for SOs, provides an additional 
problem for ACDM: such a correlation is difficult 
to explain if the formation of bulges and disks in 
SO galaxies were decoupled. 

Alternatively, SOs might be descendants of spi- 
rals whose star formation has faded after consum- 
ing the gas or losing it by some st ripping mech- 
anism , like ram press ure stripping (|Gunn fc Gott 
19721 ). halo s tripping (Sekki et al.)l2002l ). or galaxy 
harrasment (jMoore et al.l ll999). Recent evidence 
supporting this idea comes from the Tully-Fisher 
(TF) relation and from the analysis of the prop- 
erties of globular clusters in galaxies. SOs lie be- 
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low the spiral galaxie s in the TF relation , hav- 
ing lower luminosities ( Bedregal et al. 20061 ) . This 
deviation is explained by luminosity evolution of 
spiral galaxies: the transformation from spirals 
to SOs occurred at various times in the past, 
and the galaxies have been passively fading ever 
since. The globular cluster frequency (the num- 
ber of globular clusters per unit F-band lumi- 
nosity) has been used as an independent esti- 
mate of the degree to which the luminosity of SOs 
has faded relative to that of it s spiral progenito r 
( Aragon-Salamanca et al. 20061 iBarr et al.ll2007t l. 
This estimate is based on the assumption that the 
frequency of globular clusters is constant during 
the transformation process. The fact that the 
bulges in SO-galaxies also have many character- 
istics of disk-like structures, including their kine- 
matic properties (Cappellari et al. 2007), is con- 
sistent with this picture. However, if SOs were sim- 
ply passively formed from SO/a spirals it would be 
difficult to explain our finding that the fraction 
of bars is considerably lower in SOs than in SO/a 
galaxies, which are expected to be their progenitors 
in the Hubble sequence. Bars should be fairly 
robust structures, evidenced by the fact that the 
bar fraction, at least in massive luminous spirals, 
is maintained nearl y constant betwee n the redshift 
range z= - 0.84 |sheth et al.ll2008l ). 

Although the hypothesis of SOs as stripped spi- 
rals is a promising idea, an important piece of 
information is still missing in this picture. In- 
deed, bars are expected to be efficient drivers of 
galaxy evolution: the angular momentum transfer 
between gas and stars leads to gas infall and subse- 
quent star formation in the centra l regions, which 
can a dd to the mass of the bulge (jFriedli fc Bend 
1993). If the angular momentum transfer occurs 



between the bar and the halo, that leads to the 



evolution of the bar ( Debattista fc Sellwoodlioool : 
Athanassoulallioolih : a bar first grows in mass and 
length, but if the bulge mass at the same time in- 
creases due to the gas infall, that might lead to 
a subsequent weakening of the bar. Weakening of 
the bar is most efficient in strong bars with flat-to p 
surface density profiles (jAthanassoula et al.l l2005) , 
typical for early-type disk galaxies (Elmegreen & 
Elmegreen 1985). In this study we have shown 
indirect observational evidence of such evolution: 
it is tempting to think that bar weakening due to 
increased central mass concentration is the expla- 



nation for the lower fraction of bars and the larger 
fraction of ovals/lenses in SOs. We also find that 
SOs galaxies have a deficiency particularly of the 
strongest bars which fits to this picture. 

A manifestation of bar-induced secular evolu- 
tion of galaxies is probably also our finding that 
nine galaxies in our sample have massive disk-like 
bulges (< B/T > — 0.^, < n >=2.5), surrounded 
by weak bars and lenses. In principle, an increase 
in bulge mass in these galaxies could have occurred 
in a similar manner as discussed above. Also, once 
the bulge mass had increased, the bar might have 
started to weaken, leaving only a weak bar inside a 
lens. The lenses surrounding the weak bars can be 
naturally explained as relics of the evolution of the 
bar: many barred galaxies have lenses of the same 
dimension as the bar, aligned with the bar major 
axis. A problem in this scenario is that normal 
spiral galaxies don't have enough gas for making 
such massive bulges by star formation (Kormendy 
& Kennicutt 2004), not at least without a signif- 
icant accretion of extragalactic gas to the disk. 
It is still possible, however, that these galaxies 
are associated with the formation and evolution 
of bars at higher redshifts. The role of bars in 
cosmological simulations has not yet been stud- 
ied inucli J _bjjt one su ch attempt has been made 
bv lHeller etaD (120071 ). Their simulations, starting 
from initial values motivated by cosmological sim- 
ulations, include star formation, cooling and feed- 
back. In these simulations disk-like bulges form at 
early phases of galaxy evolution during the gas- 
rich epoch in the history of galaxies, being thus 
capable of accounting for the large bulge-to-total 
flux ratios of the disky bulges found in some of 
the SOs. The triaxial halos are the driving force 
to the formation of primordial bars, which trig- 
ger nuclear bars. After 4-5 Gyr the primary bars 
are weakened to fat ovals and the nuclear bars are 
decoupled. These processes might be a key for un- 
derstanding the multiple bar/oval/lens structures 
seen in many SO galaxies. 
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Fig. 1. — Two-dimensional multi-component de- 
composition is shown for NGC 524. Small white 
dots show the pixel values in the observed image 
and the other symbols show the model functions: 
dark grey lines are for the Sersic bulge (Bui) and 
the exponential disk, and small black dots show 
the final model. Lenses (L) are shown by large 
black dots and by shadowed light grey. The upper 
left panel shows the observed _ftT s -band image and 
the lower panel is the fitted model image. 
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Fig. 2.— Decomposition for NGC 5365. The 
symbols are as in Fig.l. The galaxy has two 
bars embedded in ovals, both components be- 
ing fitted by a single function (B+L). The three 
lower panels show the observed image in different 
scales. Left: the scale is selected to show the pri- 
mary bar, whereas the isophotes show the nuclear 
bar and the weak oval surrounding the primary 
bar. Middle: shows the bright inner oval, and an 
isophote indicating the nuclear bar. Right: shows 
the nuclear bar. In the ellipticity profile the two 
bar+oval systems have nearly the same ellipticity. 
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ngc3998K SA(r)0 




Fig. 3. — Left: the original if s -band image of 
NGC 3998; Right: the residual image after sub- 
tracting the bulge model obtained from our 2D 
decomposition. The bright lens at r < 12" was 
fitted by a Ferrers function. The scale is in arc- 
seconds. 
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Fig. 4. — The intrinsic scale length of the disk 
as a function of the effective radius of the bulge. 
Two if s -band magnitude bins are shown (H n =75 
km s _1 Mpc -1 is assumed). The line shows the fit 
log/i R = 0.62 + 0.46 log R eS . 
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Table 1: Bar fractions using RC3 family classes for the complete NIRSOS+OSUBGS samples. The A2 
fractions for SO and SO/a galaxies, and the statistics for ovals/lenses and multi bars have been derived from 
the NIRSOS sub-sample of 127 galaxies. 



bar-index 


SO", S0 U , S0+ 


SO/a 


Sa,Sab 


Sb,Sbc 


Sc,Scd 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


B+AB (RC3) 


46±5 % 


77±9 % 


65±7 % 


70±6 % 


61±7 % 


B (RC3) 


27±4 % 


50±10 % 


30±7 % 


39±6 % 


27±6 % 


AB (RC3) 


19±4 % 


27±9 % 


35±7 % 


31±6 % 


34±7 % 


A (RC3) 


53±5 % 


22±9 % 


35±7 % 


28±5 % 


38±7 % 


all bars (A 2 >0.1) 


53±6 % 


93±5 % 


65±7 % 


69±6 % 


64±7 % 


strong (A 2 > 0.6) 


10±3 % 


38±9 % 


26±7 % 


23±5 % 


11 ±4 % 


medium (A 2 =0.31-0.6) 


33±6 % 


44±9 % 


35±7 % 


34±6 % 


26±6 % 


weak (A 2 =0. 1-0.3) 


9±4 % 


10±5 % 


5±4 % 


11±4 % 


28±6 % 


ovals/lenses 


97±2 % 


82±9 % 


83±7 % 






multiple bars 


21±6 % 


40±12 % 


26±8 % 







(among barred) 
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